Large scale inhomogeneity of inertial particles in turbulent flow 
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Preferential concentration of inertial particles in turbulent flow is studied by high resolution 
direct numerical simulations of two-dimensional turbulence. The formation of network-like regions 
of high particle density, characterized by a length scale which depends on the Stokes number of 
inertial particles, is observed. At smaller scales, the size of empty regions appears to be distributed 
according to a scaling law. 



The transport of inertial particles in fluids displays 
properties typical of compressible motion even in incom- 
pressible flows. This is a consequence of the difference of 
density between particles and fluid. The most peculiar 
effect is the spontaneous generation of inhomogeneity out 
of an initially homogeneous distribution. The clustering 
of inertial particles has important physical applications, 
from rain generation to pollutant distribution and 
combustion Starting from the first examples in lam- 
inar flow 0, it is now demonstrated both numerically 
^ 1^ and experimentally "7, ^\ that also in fully devel- 
oped turbulence there is a tendency of inertial particles 
to form small scale clusters. The parameter characteriz- 
ing the effect of inertia is the Stokes number St, deflned 
as the ratio between the particle viscous response time 
Tg and a characteristic time of the flow Tv. In the limit 
St ^ inertial particles recover the motion of fluid par- 
ticles and no clusterization is expected. In the opposite 
limit St — » oo particles become less and less influenced 
by the velocity field. The most interesting situation is 
observed for intermediate values of St where strong clus- 
terization is observed |^ Q ■ 

In the case of a smooth velocity field the Eulerian char- 
acteristic time Tv is a well defined quantity as it can be 
identified with the inverse Lyapunov exponent — 
of fluid trajectories. In this case some general theoreti- 
cal predictions are possible 0, 0| such as the exponen- 
tial growth of high order concentration moments. De- 
tailed numerical simulations in a chaotic random flow 
have shown maximal clusterization (measured in terms 
of the dimension of the Lagrangian attractor) for a value 

St ~ 0.1 113. 

In the case of turbulent flow, where the velocity field 
is not smooth, a simple scaling argument suggests that 
maximal compressibility effects are produced by the 
smallest, dissipative scales jl^. Nevertheless, for suffi- 
ciently large values of St, the particle response time in- 
troduces a characteristic scale in the inertial range which 
breaks the scale invariance of the velocity field and pro- 
duces, as we will see, large scales inhomogeneity in par- 
ticle distribution. 

The motion of a spherical particle in an incompressible 



flow, when the size of the particle is so small that the 
surrounding flow can be approximated by a Stokes flow, 
is governed by the set of equations [l^ 
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where v represents the Lagrangian velocity of the parti- 
cle, l3 — 3po/{po + 2p) and p and po are the density of par- 
ticle and fiuid respectively. In u{x,t) represents the 
velocity field whose evolution is given by Navier-Stokes 
equations 
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In what follows, we will consider the limit heavy particles 
such that /? ~ 0. In this limit it is easy to show that the 
Lagrangian velocity possesses a compressible part |lOj : 
expanding to first order in Tg and using V • m = 0, 
one obtains, from (0) 



V • V 



sV • (m • Vu) ^ 
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From Q it is possible to give a dimensional estima- 
tion of the relative importance of the compressible part 
for a turbulent velocity field with scaling exponent h, 
Siu U{£/L)^, L and U being a characteristic large scale 
and velocity. The scaling exponent for the compressible 
component of v is 6iV ^ U{£/L)'^'^^^ and thus the rela- 
tive compressibility scales as (£/L)^~^, i.e. reaches the 
maximum value at the viscous scale [ly|. Nevertheless, 
we will see that the presence of inertial range of scales in 
the turbulent flow generates large scale structures in the 
particle distribution at large St. 

In this letter we address the problem of trans- 
port of heavy inertial particles in fully developed two- 
dimensional turbulence in the inverse cascade regime. 
High resolution direct numerical simulations of two- 
dimensional Navier-Stokes equation with white in time, 
random forcing / localized at small scales £f have been 
performed. As customary, a friction term —au is added 
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FIG. 1: Snapshots of particle concentrations taken at the same time in stationary condition for two realizations with different 
Stokes numbers St — 0.12 (left) and St = 1.2 (right) started from identical initial conditions and advected by the same turbulent 
flow. The number of particles in both cases is 1024^. Gray levels correspond to particle concentration. Turbulent velocity is 
obtained by integration of two-dimensional Navier-Stokes equation by a standard pseudo-spectral code on doubly periodic unit 
square domain at resolution A*'^ = 1024^. 



to Q in order to extract energy from the system at the 
friction scale £fi- ^ e^^^a^^^^. The intermediate scales 
<^ ( ■€i £{r define the inertial range in which Kol- 
mogorov scaling SiU U{i/Ly^^ is observed La- 
grangian tracers are placed at random with initial zero 
velocity and integrated according to with a given Tg. 
After a scratch run long several Tg, Lagrangian statistics 
is accumulated for typically some tens of Tg. Stokes time 
is made dimensionless by rescaling with the Lagrangian 
Lyapunov exponent of fluid particles, St = AiTg. Fig. 
shows typical distributions of inertial tracers in station- 
ary conditions at different values of St, obtained starting 
from the same initial random distribution. One observes 
in both cases strong inhomogeneous distributions with 
empty "holes" , in the second case on much larger scales. 

Maximum compressibility effects are expected at small 
scales and can be described by the Lyapunov spectrum 
for inertial particles. We recall that for a generic dynam- 
ical system the sum of the Lyapunov exponents gives 
the exponential rate of expansion (or contraction) of the 
hypervolume in phase space. In our case, from ^ we 
have X]i=i 4 '^j = ~2/ts, thus volumes are contracted at 
a constant rate. Let us observe that when Tg — > oo, phase 
space contraction rate vanishes, and we thus expect less 
clusterization. As a consequence of the structure of 
we find that two Lyapunov exponents are close to — 1/ts, 
representing the rate of adjustment of Lagrangian veloc- 
ity to the Eulerian one. The first Lyapunov exponent is 
found positive, as the trajectories are chaotic and the sec- 
ond, negative, determines the dimension of the attractor 
according to the definition of Lyapunov dimension |14| 



where K is defined as the largest integer such that 
^iLi ^ 0- In the inset of Fig. [5] we show the depen- 
dence of the first Lyapunov exponent on Stokes number. 
As St increases, Ai(St) decreases monotonically from the 
neutral value Ai(0) ~ 0.72. The behavior of Lyapunov 
dimension is also shown in Fig. [5] In the limit St — > 0, 
particles become neutral and thus one recovers the ho- 
mogeneous distribution with c^l = 2. The presence of a 
minimum around St ~ 0.1 was already discussed in the 
case of smooth flow and indicates a value for which 
compressibility effects are maximum. We remark that 
the curve dL{St) of Fig.|21is almost identical to the one 
obtained in synthetic smooth flow 11], and is thus insen- 
sitive to the presence of the hierarchy of scales typical of 
a turbulent flow. 

In the turbulent scenario, instead of becoming more 
homogeneous, for increasing St the inertial particle dis- 
tribution develops structures on larger scales, as is evi- 
dent in Fig. n The most evident dishomogeneities are 
related to the presence of empty regions ( "holes" ) on dif- 
ferent scales. We have thus studied the hole statistics at 
varying Stokes number. 

We have performed a coarse graining of the system by 
dividing it into small boxes forming the sites of a square 
lattice and counting the number of particles contained 
in each small box. From this coarse grained density we 
have computed the probability density function of holes, 
defined as connected regions of empty boxes. 

Probability density functions (pdf) of hole areas are 
shown in Fig. O for different Stokes numbers. The 
hole distributions follow a power law with an exponent 
— 1.8 ± 0.2 up to an exponential cutoff at a scale Ast 
which moves to larger sizes with St, as shown in the in- 
set of Fig. 131 At variance with the smooth flow case, in 
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pdf at larger scales. 

These results support the following physical picture. 
According to heavy particles filter out velocity fluc- 
tuations on timescales shorter than ts ■ We may thus ex- 
pect that inhomogeneity does appear at a scale £s whose 
characteristic time is of order of ts- For a velocity field 
with a roughness exponent h < 1, this would lead to 
Ast ^ St^^'-^^''-', i.e. to the increase of hole areas with 
St, in qualitative agreement with Fig. O Clusterization 
of inertial particles in turbulence appears to be a self- 
similar process up the cutoff scale Ast- Indeed the pdfs 
of hole areas collapse when areas are rescaled with the 
cutoff area, as shown in Fig. 21 



FIG. 2: Lyapunov dimension for heavy particles in two- 
dimensional turbulence as a function of the Stokes number. 
Inset: the first Lyapunov exponent Ai. 
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FIG. 3: Probability density functions (pdf) of hole areas, nor- 
malized with the area of the box, for St = 0.12 St — 0.6 
(x), St = 1.2 (*) and St = 2.4 (□). Holes are defined as con- 
nected regions of coarse grained distribution with zero density. 
Probability density functions are computed over 100 indepen- 
dent realizations of 1024^ particles each. In the inset we show 
the dependence of the cutoff area Ast (defined by the con- 
dition that 1% of holes has larger area) on St. The robust- 
ness of the hole area census with respect to particle statistics 
have been checked by increasing the number of tracers up to 
4 X 10^ 



which particle density recovers homogeneity for St > 0.1 
^1] , in the turbulent case inhomogeneities are pushed to 
larger scales when increasing St. 

The hole area pdf is independent on the number of 
particles used in the simulations. This is a non trivial 
property, reflecting the fact that inertial particles cluster 
on network-like structures where a clear-cut distinction 
between empty regions and particle-rich regions can be 
observed. We have also verified that the choice of the 
small coarse graining scale does not modify the hole area 
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FIG. 4; Probability density functions of hole areas rescaled 
with the cutoff areas Ast computed from Fig|3 Symbols as 
in FigEl 

A natural measure of the degree of clusterization of 
inertial particles which has been considered in the lit- 
erature 0, 0, is the deviation of the particle density 
distribution from Poisson law, D{R) — {uj^ — cr^'')/A;j, 
where A/f and cij^ represent the mean and standard devia- 
tion of particle distribution coarse grained at scale R and 
CTjj, is the corresponding standard deviation for a Pois- 
son distribution. The maximum of the deviation D{R), 
which depends on St (see Fig. 01 is then used for defining 
the typical scale of clusters, Rst- We have found that 
this definition of D{R) is strongly sensitive to particle 
statistics, nevertheless, the increase of Rst on St (inset 
of Fig. Isj seems to be a robust result. 

The presence of structures in the inertial particle dis- 
tribution is often attributed ^ to the fact that heavy par- 
ticles are expelled from vortical regions. Although struc- 
tures are related to the presence of many active scales 
in the turbulent flow, one should recall that in 2D tur- 
bulence, as a consequence of the direct vorticity cascade, 
vorticity is concentrated at the small scales [l^ and no 
large scale coherent structures appear. Holes emerge as 
a result of the delayed dynamics which filters the 
scales of the underlying turbulent flow characterized by 
times of the order of the Stokes time Tg. 
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FIG. 5: Deviation of the inertia! particle distribution from 
Poisson law as a function of the coarse grain scale for different 
St. The maximum of D{R) defines the clusterization scale 
Rst, whose dependence on St is shown in the inset. 



As discussed before, clustering occurs also in synthetic 
flows where a hierarchy of time scales is absent, just as 



a consequence of the dissipative character of the motion 
However, it appears from our simulations that to 
fully understand the geometry of inertial particle distri- 
bution in a turbulent flow the presence of structures char- 
acterized by a large set of time scales cannot be ignored. 

We conclude that the geometry of inertial particle clus- 
ters in developed turbulence is controlled both by the 
dissipative effective dynamics of the particle motion at 
small scales, and by the tendency of inertial particles to 
filter the active scales characterized by times of the order 
of the characteristic relaxation time of the particles. A 
full understanding of the geometry of particle clusters in 
developed turbulence is particularly relevant for several 
applications, such as coalescence processes or chemical 
reactions. The reaction rate of two chemical species is a 
function of their concentration. When particles of differ- 
ent species are transported by the same turbulent flow, 
their local concentrations are not independent and the 
presence of large scale correlations in the particle distri- 
bution can in principle influence the reaction velocity. 
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